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1. Introduction 


Graphical methods to study the behaviour of systems showing different chemical 
equilibrium are known and very used in several fields of Chemistry (particularly in 
Bioinorganic, Medicinal and Pharmaceutical Chemistry) in order to establish the chemical 
species that are related with drugs behaviour in different systems. Among these methods, 
the most commonly used are the distribution diagrams (Högfeldt, 1979), the titration curves 
(Asuero & Michałowski, 2011) and the molar ratio and continuous variations methods 
(Hartley et al., 1980). 


In the present work we have selected two molecules used extensively like drugs. In order to 
exemplify some novelties related with distribution diagrams and titration curves for acid- 
base systems we have selected the case of oxine (HOX, 8-hydroxiquinoline) that has been 
used as antiseptic and disinfectant. On the other hand, we have selected the complexation 
interaction between Fe(III) and tenoxicam (Tenox) to show other novelties related with more 
complicated distribution diagrams and molar ratio and continuous variations methods, 
because tenoxicam has been extensively used as non-steroidal anti-inflammatory drug that 
may be complexed with several metal ions. The chemical developed formulae of these 
compounds are presented in Scheme 1. 


2. Distribution diagrams for acid-base and complexation systems 


Graphic representations of chemical systems have found wide application because a simple 
look at them allows for solve specific problems and have a panorama, qualitative and 
quantitative, for different problems and phenomena. Moreover, some of these 
representations also permit to graphically solve the stated problem with a predetermined 
error (Vicente-Pérez, 1985). Distribution diagrams of species are some of the most used 
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graphic representations since the second half of 20th century; nevertheless, there have been 
some novelties on the field in the last decade (e.g. Moya-Hernández et al., 2002a, 2002b). 


OH O LK 


OH 
(a) (b) 


Scheme 1. Developed formulae of a) oxine (HOX) and b) tenoxicam (Tenox). 


2.1 Typical representations of distribution diagrams 


In this section some typical representations will be quickly reviewed to introduce the newest 
representations given in section 2.2. 


2.1.1 Acid-base systems 


In Chemistry the study of polyprotic systems, whose global formation equilibrium are 
represented in Eq. 1, is of crucial importance, because many substances, drugs among them, 
follow this Bronsted acid-base behaviour. 


, (LH) *] 
L® + jH* => LHI” with g,=——1__ (1) 
i 11H y) 
The typical way to define the molar fractions to describe the species distribution of the 
component L in this system with respect to the proton (H+) is given in Eq. 2, as well as its 
factorization in the substance amount balance equation of this component in the system - 
using the set of Eq. 1 (Rojas-Hernández, 1995). 


LAY AY 


j=1 


f= (2) 


where [L]ris L total concentration in the system. 


An example of distribution diagram is given in Fig. la, showing the case of oxine 
hydrochloride (H2OXCI, 8-hydroxiquinolinol = HOX). 


Fig. 1b represents the function known as Average Proton Number (7), introduced and 
developed by Niels and Jannik Bjerrum (Hartley et al., 1980). Eq. 3 shows the definition and 
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factorization of n as a function of pH (taking into account the sets of Eq. 1 and 2) for 
aqueous solutions. 


y jA IHY 
py 2 n 


n= = = 


i if; (3) 
[L]; 1+5 [H Ý j=1 á 
j=1 


[H*]r is the protons total concentration in the system, which requires the balance of proton 
equation. In aqueous solutions, this equation needs the autoprotolysis constant and special 
considerations. 


The chemical information given in the curves of Fig. 1 is the same, as expected from Eq. 2 
and 3. In fact, the Average Proton Number is the set of statistical means of the subjacent 
distributions resumed by the distribution diagram of Fig. 1, as it will be explained in the 
section 2.2. 
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Fig. 1. Typical graphic representations of the oxine hydrochloride (H20XCI) in water as a 
function of pH. j (or v) represents the proton stoichiometric coefficient for each species. 

a) Distribution diagram. b) Average proton number. Data were taken from Ringbom (1963): 
pKa1 =5.0, pKaz = 9,7. 


2.1.2 Complexation systems 


Even though the knowledge of Bronsted acid-base behaviour is important, the interaction of 
substances with metal ions is remarkable as well (e.g., the interaction of drugs with metal 
ions may potentiate or suppress its pharmacological activity or toxicity). 


The formation of several coordination compounds, or complexes, between a metal ion (M) 
and a ligand (L) can be described by the global formation equilibrium represented in Eq. 4. 


bos _ My 
iM+jL a NAPS with Bij Sn (4) 


where ie (1, 2, ..., m} and je {(0, 1, ..., n[i]) 
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In Eq. 4, the charge of the species has been omitted for notation simplicity. When i = 1 the 
complexes are called mononuclear, but when i 2 2 the corresponding complexes are called 
polynuclear. The formation of polynuclear complexes on a given system is 
thermodynamically favoured when the total concentration of M ([M]r) is high and placed 
over the mononuclear wall (Ringbom, 1963). 


In complexation chemistry, several distributions have to be considered. In general, it is 
preferred to study the way the component L is distributing on M species, but studying how 
the component M is distributing on L species may be interesting as well. In both cases there 
are two possible descriptions, depending on what is considered between the amount or the 
concentrations of the species. When polynuclear compexes are formed in the system all the 
distribution diagrams representing the formation of the species depends on [M]r and [L]r. 


Distributions of L in M 


The M substance amount fractions are defined by Eq. 5. 


t/M¢ — [M] _ 1 
fio [Mr | Dp mey | 
i=1\j „Pil (5) 
[M;L;] ra 
and Mf, = a = /M GRIMIT) 


where ie {1, 2, ...,m} and je {0, 1, ..., n[i]} 


When polynuclear M species are not forming in a given system (i = 1), this distribution only 
depends on [L]r, but when polynuclear species appear, the substance amount fractions 
depend on [L]r and [M]7; furthermore in this last case the simple sum of the concentrations 
of M species (Em) is lower than [M]r, in agreement with Eq. 6. 


m { ni] ni] 
Em = $m 11) Zama, 1|- [MJ (6) 


i=1\ j=0 i=1\ j=0 


In this particular case, M concentration fractions can be defined by Eq. 7. 


itty, 0) - Mg E Mop IMPLY) O) 


= L and Q; = 


where ie {1, 2,..., m} and je (0, 1, ..., n[i]} 


The substance amount and concentration fractions of M species are related by means of 
Eq. 8. 


L/M j e eN (8) 
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Distributions of M in L 


Following the same approach, several distributions of L species may be defined. Then, the L 
substance amount fractions are defined in Eq. 9, while the L concentration fractions are 
described in Eq. 10. 


M/L p= [L] _ 1 
[L]r a (W Lye? 
14> > ¡8 [MIIL ^ 
j=l i=0 (9) 
and MI" f,= PE a po JB MITIA) 
[L]y 
where ie (0, 1, ..., m[j]} and je (1,2, ..., n} 
- [L] _ 1 
M/L Py B [ Jr ~ i 2 [M] [LJ 1) 
+ iB; 
j=1\ i=0 ý (10) 
and Mig = LOM) meo Cp IMEL”) 
[L]r 


where ie {0, 1,..., m[j]} and je {1, 2, ..., n} 


The substance amount and concentration fractions of L species are related by means of 
Eq. 11. 


mig 1672, lik an 
L 


because the following inequality is always confirmed. 


n (mljl n ( mj] 
2, = (Ptos) aps fst} =1Lh (12) 


jali 
Distribution diagrams for the Fe(III)-tenoxicam system in acetone 


In order to exemplify the typical distributions diagrams of M and L species, Fig. 2 is 
reported for the Fe(III) - tenoxicam (Tenox) system in acetone from previously reported data 
by Moya-Hernández et al. (2009). The equilibrium constants of this reference have been 
collected in Table 1. 


| Species |iļj| log% | 


14.75 + 0.06 
18.45 + 0.07 


Table 1. Global formation constants of Fe(III)-tenoxicam species in acetone (Moya- 
Hernandez et al., 2009). 
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The distribution diagrams that represent the substance amount fractions, calculated from 
global formation constants given in Table 1, have been constructed with the aid of program 
MEDUSA (Puigdomenech, 2010). The distribution diagrams that represent the concentration 
fractions have been obtained by means of Excel (Microsoft®) worksheets applying Eq. 8 and 
11. Some of the typical distribution diagrams are shown in Fig. 2 and 3. 


Fe,Tenox 

10 Fe,Tenoxs Fe 1.0 e Fe 

0.8 0.8 5 
eS S 
S o5 $ os] 
2 3 
ES & 

0.3 4 0.3 4 

0.0 

0.0 > 
1.0 2.0 3.0 4.0 5.0 
1.0 2.0 3.0 4.0 5.0 
pTenoxr pTenox; 
(a) (b) 


Fig. 2. Typical graphic representations of the Fe(III)-tenoxicam system in acetone as a 
function of pTenoxr. [Fe(11)]r = 1x103 M and pTenoxr = -log[Tenox]r. a) Distribution diagram 
of substance amount of Fe(III) species taking into account the quantity of Fe(III) in each 
species. b) Distribution diagram of concentration of Fe(III) species taking in account which 
Fe(III) species is more concentrated in the system. 
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Fig. 3. Typical graphic representations of the Fe(III)-tenoxicam system in acetone as a 
function of pTenoxr. [Fe(III)]r = 1x103 M and pTenoxr = -log[Tenox]r. a) Distribution diagram 
of substance amount of Tenoxicam species, taking in account the quantity of tenoxicam in 
each species. b) Distribution diagram of concentration of tenoxicam species, taking in 
account which tenoxicam species is more concentrated in the system. 
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The independent variable of these figures, in all cases, has been selected to be the pTenoxr, in 
order to compare directly the shape of each fraction in Fig. 2 and 3. Nevertheless this axis 
could represent another variable (like pX, p21, pL, pM, etc.) 


It is noteworthy than distribution of substance amount and concentration are similar for the 
same component, even though slightly differences can be shown comparing Fig. 2a and 2b, 
or Fig. 3a or 3b; some of their differences are crucial and related with the physical-chemical 
meaning of each distribution. 


A distributions of Fe(III) and tenoxicam species comparison allows for the conclusion that 
they are very different, obviously representing distributions of two different components in 
the system. 


2.2 Distribution diagrams as a function of stoichiometric coefficients 


The sentence “distribution diagram” could have an implicit idea concerning possible 
statistical distributions subjacent to the graphic representations given in Fig. 1, 2 and 3. This 
idea was first explored by Moya-Hernández et al. (2002a) for the case of Bronsted acid-base 
systems. In the next subsection (2.2.1) the conclusions of this work will be applied to the case 
of oxine species as a function of pH while in the other subsection (2.2.2) this treatment will 
be generalized to the case of complex systems of the M-L kind, where polynuclear species 
are forming. 


2.2.1 Distribution diagrams of one discrete variable (the proton stoichiometric 
coefficient for an acid-base system species ) 


If an aqueous solution of oxine at a certain concentration has a known pH value, the 
distribution of the oxine species will be fixed, as it is shown in Fig. la. Then a 3D graph 
could represent all the distributions of oxine species for each pH value. The set of 
distributions of oxine species is represented in Fig. 4 as well as one of them, at pH = 5.0. 


In agreement with Mathematical Statistics (Kreyszig, 1970; Reichl, 1980), the meaning of 
each distribution of discrete variable is defined by Eq. 13. 


v=) if (13) 
j=0 


Furthermore, each of these statistical distributions has a variance. The variance of each 
distribution of discrete variable is defined by Eq. 14. 


=D VIS (14) 


The equality of Eq. 3 and 13 demonstrate that the set of means of the oxine species 
distributions is the one given in Fig. 1b. In this way, the mean of the distribution offers the 
value in which the proton number is centered, as an average, at each pH value. 


The set of variances of these distributions is related with an intrinsic buffer capacity (Moya- 
Hernández et al., 2002a). A graphic representation of the oxine set of variances as a function 
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of pH is given in Fig. 5. When the variance has a value near to zero one species is present in 
the system with a fraction near 1; when the variance reaches a maximum, two or more 
species are present in the system with comparable fraction values. In the case of oxine 
system the maxima of the variance function are reached for pH values equal to pK, and 


PKa2, because pKa >> pKa. 


HOX E i 1.0 
pH= pKa =5.0 


Molar fractions 
Molar fractions 


(a) (b) 


Fig. 4. Statistical distributions of oxine species. a) Set of distributions of the stoichiometric 
coefficient of protons as the discrete variable. b) Specific distribution diagram at pH = 5.0. 


0.3 


Set of variances 


Fig. 5. Set of variances of the distributions of discrete variable of oxine as a function of pH. 
The maxima of these function are placed in the pK, values. 


2.2.2 Distribution diagrams of two discrete variables (the M and L components 
stoichiometric coefficients for complexation systems) 


It can be demonstrated that the distribution diagrams as those of Fig. 2 and 3 for M-L 
complexation systems represent sets of distributions of two discrete variables,where these 
variables are the M and L stoichiometric coefficients. 
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Examples of substance amount and concentration distributions of tenoxicam in Fe(III), for 
only one of the systems represented in Fig. 2 are shown in Fig. 6. 


In this case, each one of the statistical distributions that can be defined has two means, two 
variances and one covariance (Reichl, 1980). 


Just as a mere example, the definition of the two means of the substance amount distribution 
of Tenoxicam in Fe(III) is given in Eq. 15, while the definition of the two variances and the 
covariance for the same distribution is given in Eq. 16. 


m ni] 
ua > iS 
i=1j=0 
n mj] (15) 
MT = Al 
j=0i=1 


Fe,Tenox, 


Tenox/Fe Qij 


Teo. HE Vre 


Fig. 6. Statistical distributions of tenoxicam in Fe(III) species in acetone. [Fe(111)]7 = 1x103M 
and [Tenox]r = 6x104M. a) Substance amount fractions of Fe(III) species. b) Concentration 
fractions of Fe(III) species. 


The set of means and variances for the distributions of two discrete variables of Tenoxicam 
in Fe(III), represented typically in Fig. 2a like a “distribution diagram”, are shown in Fig. 7. 


The interpretation of each of the two means is to be the stoichiometric coefficient average for 
the corresponding component. The two variances grow when two or more species are 
forming: higher variance increases the number of species with different stoichiometric 
coefficient of the corresponding component. 


The possible consequences of this statistical view should be studied exhaustively, but this 
kind of study would be beyond the objectives of the present work. The distribution 
diagrams applications to determine or to use stoichiometric coefficients are developed in the 
following sections. 
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Fig. 7. Statistical parameters for the substance amount distributions of Tenoxicam in Fe(III) 
in acetone corresponding to Fig. 2a. a) Set of the two means of the distribution. b) Set of the 
two variances of the distribution. 


4.0 5.0 


3. Acid-base titration curves for polyprotic systems 


This section deals with titration curves, pH = f(volume of strong acid or base), of polyacid 
and polybase systems by means of a description given by a thermodynamic model. This 
allows for the theoretical building of this kind of curves as well as their first derivative and 
the buffer capacity curves as a function of pH. The model is deducted from electroneutrality 
equation relating the global formation equilibrium and molar fractions of the distribution 
diagram of the system. Some previous works concerning this mater can be found in 
literature, by Fleck, 1967; Högfeldt, 1979; King & Kester, 1990; Efstathiou, 2000; Rojas- 
Hernández € Ramirez-Silva, 2002; Tarapcík & Beinrohr, 2003; Asuero, 2007; Gutz, 2010; 
Asuero € Michatowski, 2011. 


3.1 Titration curves, pH = f(added volume of strong acid or base), of a mixture of 
species in the same polyprotic system 


A system formed by a mixture of V,¡ volumes of solutions of HL“ species, with 
concentrations Coj, giving a total initial volume V, = XV., has been considered. The HjLi« 
species form part of the same polyprotic system (H,L(-W*/H (1.1 L(a-D3/.../HaL/.. /[HLG-D;/Lo- 
/ H+). Each charged species has associated a contra-cation or a contra-anion (Mz* or Xz) that 
do not have acid-base properties, depending if (j-a) is negative or positive. 
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3.1.1 Titrations with a strong base: MOH 


If the mixture described at the beginning of section 3.1 is titrated with a strong base MOH, of 
Cp concentration, it is possible to write the electroneutrality equation, for each added volume 
of MOH (V?+), in the form: 


a-1 
ASA 
+ I3 j-a + 
[M*]+ AENA + LU a)[LH}“]+[H*]= 
> (17) 
IA 
_ j=a+1 i j-a L 
= Tay, * La j[LH}“]+[OH] 


where [M+] = V,C,/(V, + Vb), the second term in the first member of Eq. 17, represents the 
contra-cations charge associated to the anions of the polyprotic system for j e (0, 1, ..., a-1}, 
while the first term in the second member represents the contra-anions charge associated to 
the cations of the polyprotic system for j e {a+1, a+2, ..., n}. 


Introducing Eq. 1 and 2 in Eq. 17 and algebraically rearranging it, it is possible to arrive to 
Eq. 18. 


j=0 


ONE c=] Emo AE g-an- 


V, = 


(18) 


3.1.2 Titrations with a strong acid: HX 


If the mixture described at the beginning of the section 3.1 is titrated with a strong acid HX, 
of C, concentration, it is possible to write the electroneutrality equation, for each added 
volume of HX (V,), in the form: 


a-1 
2 (1 PAVGCoj} 
j=0 o j-a +7 
yy Dee I 
n (19) 
> (j ~ AVC} a-1 
_ j=a+1 o; j-a = — 
= ATA * La DILH?"]+[0H7]+[X7] 


where [X] = V,C,/(V. + Va). Introducing Eq. 1 and 2 conveniently in Eq. 19 and 
algebraically rearranging it, it is possible to arrive to Eq. 20. 
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Fig. 8. Titration curve of 50 mL of aqueous solution of H:OXCI with an aqueous solution of 
NaOH 0.1315 M. The values used to obtain the fitting shown are the following: pKa = 5.0, 
pKa = 9.7, pKw = 13.7, [HMOXCI]initiar = 0.00101 M. 

1) The circles represent the experimental points of the pH = f(V,) curve and the solid line is 
the fitted curve obtained with Eq. 18. 2) The rhombuses represent the experimental points of 
the (dpH/dV;) = g(V,) curve and the dashed line is the fitted curve obtained with Eq. 22. 

a) Curves as a function of added volume of the titrand. b) Curves as a function of molar 
ratio (r) of the titrand. 
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Practically all the simulators, available nowadays, to predict acid-base titration curves, using 
strong base or acid as titrand agent, are based on Eq. 18 and 20. They are analytical 
expressions to calculate exactly the added volume of strong base or acid, given a set of pH 
values. In the following subsections this feature will be used to obtain important 
applications of these equations. 


-D-a Vj Cy) Es oj Coj | [= a) fi} jw Pa 
V, = (20) 
C,-[H* |» Ky 
[H"] 


3.1.3 Titration of aqueous solution of oxine hydrochloride (H20XCI) with NaOH 


In order to compare the predicted results with the experimental ones, 50 mL of an aqueous 
solution 0.001M H2OXCI where titrated with a solution of NaOH 0.1315 M at 25°C. 


The comparison of the experimental titration curve with the fitted curve through the model 
given by Eq. 18 is shown in Fig. 8a. 


As it is shown in curves 1 of Fig. 8, the observed fitting for the titration curve is good. In the 
case of Fig. 8b, the molar ratio (r) is defined as the ratio of the added titrand (NaOH), with 
respect to the analyte (OX7r). The solution pH has a great change near 1, meaning that the 
quantitative reaction consumes 1 mol of H20X* for each added mol of OH”. In other words, 
the stoichiometric coefficients of these reagents are 1 and 1 for the reaction that permits the 
quantification of oxine with the hydroxide ion: 


H,OX* +OH ——HOX+H,0 (21) 


3.2 The first derivative of the titration curve: dpH/dV 


The first derivative method of an acid-base titration curve is well known to determine the 
equivalence point volumes in order to accomplish quantitative chemical analysis. 


Nevertheless, being Eq. 18 and 20 functions of only one variable (the pH), it is possible to 
obtain the analytical expressions of their first derivatives and their reciprocal functions, to 
arrive to the exact algebraic expressions of the titration first derivatives curves. These 
expressions are given in Eq. 22 and 23. 


-pH = pH-pK, 
dpH _ C, +10 10 22) 
dv, Z z H H-pK 
-2.303 De Co) NOIA] | + 2:303(V, + V,)| 1077" +107 | 
jol is 
= -pH pH-pK, 
dp H _ C,=107" +10 (23) 
dV n 
-2.303| X (VC) || > OS +2.303(V, +V,)[ 107" + 1075. | 
j=0 ¡=0U  i=0 
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3.2.1 Titration first derivative curve of an oxine hydrochloride (H20XC/) aqueous 
solution with NaOH 


The curves 2 in Fig. 8 compare the titration experimental first derivative, described at the 
beginning of subsection 3.1.3, with the curves obtained by Eq. 22. 


The experimental first derivative was obtained approximately as the ratio of finite 
differences of measured pH values and volumes during titration, using the average of 
volumes, or molar ratios, for each interval. 


As it can be seen, the fitting attained is quite good and the maximum observed in Fig. 8a and 
8b is sharp. In the first case, this maximum signals the volume position of the first 
equivalence point, while in the second case it indicates the ratio of the titrand stoichiometric 
coefficients with respect to the analyte for the quantitative reaction. 


The first derivative is usually better than the pH curve to experimentally determine the 
volumes of the equivalence points (Fig. 8a) and the titration reactions (titrand/analyte) ratio 
of stoichometric coefficients (Fig. 8b) when these are quantitative. 


In the case of the HOX reaction with OH, the second titration reaction, is not quantitative 
because its corresponding equilibrium constant is not high enough for a 0.001 M analyte 
initial concentration. For this reason there are no visible changes in the second equivalence 
point volume in curves 1 or 2 in Fig. 8a, nor in Fig. 8b for r = 2. 


3.3 The buffer capacity (8) of a polyprotic system 


In many chemical and biological processes it is essential that the medium pH be kept within 
certain limits, which is possible through the use of buffer solutions. They possess a specific 
buffer capacity and are used to maintain constant the pH with a very small uncertainty. 


In the chemical literature, there are two ways to define a buffer capacity (PB): one is defined 
in terms of concentration of strong acid or base added to the system, in order to simplify the 
mathematical treatment, as firstly proposed by Van Slyke (1922), and then used by others, as 
Urbansky & Schock (2000) or Segurado (2003). The other way to define the buffer capacity is 
in terms of the amount of strong acid or base added to the system, as King and Kester did 
(1990), as well as Skoog et al. (2005); they also derive equations with the concentration, but 
considering explicitly 1L of solution. 


3.3.1 A buffer capacity with dilution effect (Bain) 


According to the definition given by King & Kester (1990) it is possible to apply the 
derivative of the added amount of strong base or acid with respect to the pH, and then 
obtaining mathematical expressions for the buffer capacity as function of pH, by considering 
the dilution effect, i.e. Bai = f(pH). In the present work, this implies to take Eq. 22 and 23 
reciprocals and multiply them by C; or C,, respectively. 


The typical way to represent a buffer capacity curve, consists in plotting it as a function of 
pH, even though it may be represented as a function of the titrand volume or the molar ratio 
titrand/ analyte. 
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Baa. = AVCy _ 2.303C, 
dil = T 
b è dpH  C,+(10”" -10 PKw) led 
[Esc E (Eso 2) +(V, +Vp)| 107" 410° -PKw | 
j=0 j=0 i=0 
Ban. = dV,Ca _ 2.303C, 
dil = — 
a dpH Es = (10 pH _ 10H PK) 05) 


j=0\  i=0 


[gr] Sendo. wm 


As Moya-Hernández et al. (2002b) have demonstrated, Eq. 14 is equal to the double sum 
between brackets in Eq. 24 and 25. For this reason, the set of variances is related with the 
intrinsic buffer capacity of a polyprotic system. 


n n n 
= 6-0) f= Sny - asi) = intrinsic buffer capacity (26) 
j=0 j=0\ — i=0 

It is noteworthy that the intrinsic buffer capacity, i.e. the set of variances of the proton 
stoichiometric coefficient distributions for this kind of systems, is a function that only 
depends on the molar fractions and the stoichiometric coefficients. 


The second term between keys in Eq. 24 and 25 is due to the acid and basic particles of the 


amphiprotic solvent (H+ and OH” in the case of water). (See Segurado (2003), Urbansky & 
Schock (2000).) 


3.3.2 Buffer capacity with dilution of aqueous solution of oxine hydrochloride 
(H20XCI) with NaOH 


Fig. 9 shows the comparison of experimental buffer capacity, with effect of dilution, and the 
curve obtained by Eq. 24 and 25 for the titrations of the system defined at the beginning of 
subsection 3.1.3. 


Bay /mmol 


Fig. 9. Comparison of experimental (markers) and calculated (solid line) buffer capacity with 
effect of dilution. 
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The experimental Bai has been determined approximately by taking the ratio of the titrand 
concentration and the first derivative experimentally obtained as explained in subsection 
3.2.1. 


As it can be seen in Fig. 9, there is a good agreement between the theoretical and the 
experimental data. The expected maximum of the curve at pH = pKx = 9.7 (see Fig. 5) is lost 
due to the effect of the hydroxide ion over the pair HOX/OX” by the low concentration of 
oxine in the system (0.001 M). 


4. Molar ratio and continuous variations methods for complexation systems 


For complexation reactions the two traditional methods to determine the stoichiometry of 
formation reactions and formation equilibrium constants are the molar ratio and the 
continuous variations methods (this last also known as Job’s method). (See Hartley et al., 
1980.) 


Nevertheless, the explanation of these methods has not always been given in a clear way, 
and this is particularly evident for equilibrium considering the formation of polynuclear 
species. Furthermore, these treatments are disappearing from modern text books. For this 
reason, the aim of this part of the chapter is to describe a way to deduct the equations and 
the curves that appear in some Inorganic and Analytical Chemistry books. 


4.1 Tables of Variation of Substance Amount (TVSA) 


Professor Gaston Charlot (1969), in France, has used a teaching tool to explain the advance 
degree in the formation of acid-base, complexation and redox reactions. It has been called by 
him the Table of Variation of Substance Amount. In the following subsections this tool will 
be applied to the systems equilibrium states using the molar ratio and continuous variations 
methods. 


For a system where only one reaction takes place, quantitatively: 


[M,Ly] 


[M]}" [L]" 


>>1 (27) 
4.1.1 Table of Variation of Substance Amount for the molar ratio method 


In the molar ratio method several amounts of one of the reagents, i. e. L, is added to another 
reagent, i. e. M, and typically the dilution of the solutions is controlled in volumetric flasks 
of the same capacity Vr in order to measure a response R that is directly proportional to the 
reaction product, M,L», through the equation: 


R=ky 1, [M,L,]= kun, ("M,L / Vr) (28) 


This experiment is described in Table 2 for a quantitative reaction. 


If all the amounts in the first column of Table 2 are divided by the M quantity (nm), and the 
amounts in the other columns are divided by Vr and multiplied by the corresponding 
response factor, Table 3 will be obtained. 
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aM + bL M,Lp 
Initial NM 
b 

ny, < Mu nm- (a/0)n;, (1/b)n1 
b 
—Ny = sc" ae =0 (1/a)nm 
a 

b 

qa <n, ni - (b/a)nm (1/a)nm 


Table 2. TVSA for molar ratio method. ny is a constant for all the systems, nz is variable, all 
systems are in solution in thermodynamic equilibrium and the total volume is the same and 
equal to Vr. scn; represents the amount of L for the stoichiometric condition. 


The functions deducted for molar ratio method and presented in Table 3, with the given 
restrictions, lead to represent their behaviour in Fig. 10 assuming Eq. 28 and giving R in 
arbitrary units. 


aM + bL Ss MaLo 
Initial NM 
b 
Be km(nm- (a/b)n1)/Vr KmaLo(1/0)r (nm)/Vr 
b 
sch == ae =0 kmare(1/a)nm/ Vr 
a ar ae =0 ki (nz - (b/a)nw)/Vr kmare(1/a)nm/ Vr 


Table 3. TVSA for molar ratio method expressed in terms of molar ratio of L (rz) and 


response (R) of the system. scr; represents the amount of L for the stoichiometric condition. 
For this study it has been assumed that ky = k; = 0. 


It should be noticed in Fig. 10 that the intersection between the two formed straight lines 
indicates the molar ratio corresponding to stoichiometric conditions; in other words, this 


value is equal to the L stoichiometric coefficient ratio divided by the M stoichiometric 
coefficient : scr, = b/a. 


(a) 


Fig. 10. Curves Response = R= f(r) for experiments of molar ratio method. a) a = 1 for the 
formation of species of the ML; kind. b) a = 2 for the formation of species of the M2Lp kind. 
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If some of the assumptions formulated in this work are avoided there will be some 
deviations of the curves shown in Fig. 10. As an example, if more species linearly contribute 
to the response, the parameters of the straight lines change, but their intersection still 
indicates the molar ratio of the stiochiometric conditions. If the formation reaction in Eq. 27 
is not quantitative there will be some deviations of the linearity near the stoichiometric 
conditions. Finally, if there are more species forming in the system at the same time, in 
principle, there will be more straight lines in the molar ratio curve, one for each species 
appearing in the system (if all the reactions are quantitative, if all the contributions to the 
response are linear and the response factors are different enough). 


4.1.2 Tables of Variation of Substance Amount for the continuous variations method 


In the continuous variations method several amounts of the reagents, L and M, are mixed in 
such a way that the total amount of both reagents is constant: nm + n, = nr = constant. 
Typically, the solutions dilution is controlled in volumetric flasks of the same capacity Vr in 
order to measure a response R that is directly proportional to the reaction product, M,Lp, 
through Eq. 28. This experiment is described in Table 4 for a quantitative reaction. 


Generally, in this method the equilibrium conditions are expressed in terms of the molar 
fraction of one of the components, x; or xy. So the scx (the component L molar fraction at 
stoichiometric conditions) can be deducted, as it is shown in Eq. 29. 


b n b 
= lt > (nr = sen) = scn, 
b 
b Z 
> sex, =—(1- sc X1) => sex, = a b (29) 
a 
1+- 
a 
x b 
SCL q+b 
nm +t ny = NT aM bL MaLo 
Initial NT - NL nL 
Ny < sc nr- ni - (a/b)mi be =0 (1/bJni 
ni Pet ae =0 be =0 (1/b)( senı) 
sc”, <n, ae =0 ni - (b/a)( nr- nı) (1/a)( nr -nz) 


Table 4. TVSA for continuous variations method. nr is a constant for all the systems, nm and 
nı are variable, all systems are in solution in thermodynamic equilibrium and the total 
volume is the same and equal to Vr. scn, represents the amount of substance of L for the 
stoichiometric condition. 


If all the amounts in the first column of Table 4 are divided by nr, while the amounts in the 
other columns are divided by Vr and multiplied by the corresponding response factor, Table 
5 will be obtained. 


The functions deducted for continuous variations method and presented in Table 5, with the 
given restrictions, lead to represent their behaviour in Fig. 11 assuming Eq. 28 and giving R 
in arbitrary units. 
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xm+xL=1 aM + bL 5 M,Lp 
Initial Nr - Ny ny 
b 
(<a oa km(n7-((b-a)/b)n1)/Vr be =0 (Kmaw0)xin1/Vr 
b 
Wa = sc* ae =0 be =0 (kmatt/(atb))n1/Vr 
a+b 
b 
ae <x <1 ae =0 kı(nı-(b/a)(nr-n1)/Vr (kmatw/a)(1-x1)n1/Vr 


Table 5. TVSA for continuous variations method expressed in terms of L molar fraction (xr) 
and response (R) of the system. scx, represents the amount of L substance for the 
stoichiometric condition. For this study it has been assumed that km = k; = 0. 


It should be remarked in Fig. 11 that the intersection between the two formed straight lines 
indicates the molar fraction corresponding to stoichiometric conditions; in other words, this 
molar fraction is equal to the following relationship: scx, = b/(a+b). 


Moreover, in this method, as it has been discussed for the molar ratio method, if some of the 
assumptions formulated in this work are avoided there will be some deviations of the 
curves shown in Fig. 11. As an example, if more species contribute to the response, the 
parameters of the straight line change, but their intersection still indicates the molar fraction 
of the stiochiometric conditions. If the formation reaction in Eq. 27 is not quantitative, there 
will be some deviations of the linearity near the stoichiometric conditions. Finally, if there 
are more species forming in the system at the same time, in principle, there will be more 
straight lines in the molar ratio curve, one for each species appearing in the system (if all the 
reactions are quantitative, if all the contributions to the response are linear and the response 
factors are different enough). 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
XL XL 
(a) (b) 


Fig. 11. Curves Response = R= f(x) for experiments of continuous variations method. a) a = 1 
for the formation of species of the ML; kind. b) a = 2 for the formation of species of the M2Lp 
kind. 
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4.1.3 Some remarks for the use of Tables of Variation of Substance Amount 


The TVSA are a very powerful tool to understand the composition of thermodynamic 
equilibrium states when chemical reactions are involved. In the subsections 4.1.1 and 4.1.2 
the dependent and independent variables are those typically selected. Nevertheless other 
selections could be used if the description and interpretation of chemical behaviour require 
it. Then, in the molar ratio method the representation of the curves could be done as a 
function of metal ion ratio (rm), while for the continuous variations method the molar 
fraction of metal ion (xm) could be used. 


Furthermore, the TVSA could be used for the treatment of several complexes forming in the 
system at the same time. 


In this work we have presented only the curves that may help interpreting the system that 
will be treated in the following section. 


4.2 Application of the methods to determine the stoichiometry and equilibrium 
constants: The case of Fe(Ill)-tenoxicam system in acetone 


A spectroscopic study in the visible region of electromagnetic spectrum was undertaken to 
determine the complexes formed as well as their formation constants for the Fe(III)- 
tenoxicam system in acetone, due to the red colour observed when both reagents are mixed. 
The solutions preparation and the details of the spectra acquisition have been given 
elsewhere (Moya-Hernández et al., 2009). 


The absorption spectra obtained in this study are presented in Fig. 12, for molar ratio and 
continuous variations methods. 


In order to be near the fulfilment of the conditions selected in section 4.1, (e.g. absorption 
only due to the complexes formed) a wavelength of 520 nm was selected to represent the 
typical curves of molar ratio and continuous variations, which are shown in Figure 13. 
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Fig. 12. Absorption spectra in the visible region for the Fe(III)-tenoxicam system in acetone. 
[Fe(III] = 0.001M or [Fe(III] + [Tenox] = 0.001M. a) Molar ratio method. b) Continuous 
variations method. 
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Fig. 13. Spectral behaviour, at 520 nm, of the Fe(III)-tenoxicam system in acetone where 
different complexes are forming. [Fe(III)] = 0.001M or [Fe(III)] + [Tenox] = 0.001M. a) Molar 
ratio method. b) Continuous variations method. The dashed lines and the thin vertical lines 
only have been added to help follow the discussion given in the text. 


The shape of the curves presented in Fig. 13 demonstrates the formation of polynuclear 
species in the system, both methods clearly showing an abrupt change of slope before the 
molar ratio or molar fraction corresponding to 1:1 stoichiometry (rtenox < 1 and 
XTenox < 0.5), as it can be concluded form Fig. 10b and 11b. The other changes in slope 
support the possibility of the formation of several complexes; in other words, two straight 
lines are not enough to explain the shape of the curves. The thin vertical lines indicate the 
position of the molar ratio or the molar fraction that fulfil the stoichiometric conditions of 
the Fe(III)-tenoxicam complexes (FezTenox, FezTenox and FezTenox3) obtained in the 
treatment of spectrophotometric data (Moya-Hernandez et al., 2009 and Table 1). 


In order to better describe the physical-chemical behaviour observed in the molar ratio 
method, the distribution diagrams of discrete variable of tenoxicam in Fe(III) (for the 
concentration of Fe(III) species) for the thin vertical lines shown in Fig. 13a, are presented in 
Fig. 14. 


As it can be seen in Fig. 14, the predominant complex in the system corresponds to the 
molar ratio of each system, with other Fe(III) species been important as well. This is the 
reason why the experimental points are placed outside straight (dashed) lines in Fig. 13a. 


Although some of these distribution diagrams of discrete variable could be represented for 
the continuous variations method, Fig. 13b may be the best graphic representation could be 
a distribution that consider Fe(III) and tenoxicam species at the same time, being the sum of 
substance amount of Fe(III) and tenoxicam (nfe +N Tenox = NT = constant) the biggest restriction 
of this method. A look of this kind of distributions have been previously developed (Moya- 
Hernandez, 2003). 
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Fig. 14. Distribution diagrams of concentration of tenoxicam in Fe(III), for Fe(III) species, of 
some of the systems related with Fig. 13a. a) fTenox = 0.5. b) rtenox = 1.0. C) rrenox = 1.5. 


5. Conclusion 


In this work, some novelties related with the distribution diagrams, understood as statistical 
distributions of discrete variables, have been presented. We have extended this vision for 
two component systems. Moreover, the fractions of the distribution diagrams have been 
used to develop analytical equations to calculate exactly pH-metric titration curves as well 
as their first derivatives and a buffer capacity with effect of dilution. Finally the molar ratio 
and continuous variations methods have been reviewed, according to Charlot’s methods, 
and explained with the discrete variables statistical distributions for complexation systems. 
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